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ABSTRACT 

Accurate estimation of cosmological parameters from microwave background aniso- 
tropics requires high-accuracy understanding of the cosmological model. Normally, a 
power-law spectrum of density perturbations is assumed, in which case the spectral 
index n can be measured to around ±0.004 using microwave anisotropy satellites such 
as MAP and Planck. However, inflationary models generically predict that the spec- 
tral index n of the density perturbation spectrum will be scale-dependent. We carry 
out a detailed investigation of the measurability of this scale dependence by Planck, 
including the influence of polarization on the parameter estimation. We also estimate 
the increase in the uncertainty in all other parameters if the scale dependence has 
to be included. This increase applies even if the scale dependence is too small to be 
measured unless it is assumed absent, but is shown to be a small effect. We study the 
implications for inflation models, beginning with a brief examination of the generic 
slow-roll inflation situation, and then move to a detailed examination of a recently- 
devised hybrid inflation model for which the scale dependence of n may be observable. 
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1 INTRODUCTION 

One of the most important legacies of the COBE satellite 
for the inflationary cosmology was its emphasis that observa- 
tions had reached a quality where the inflationary prediction 
for the density perturbations could no longer be taken to 
be the scale-invariant Harrison-Zel'dovich spectrum. That 
paradigm was quickly replaced by a new one, where the 
spectrum is approximated by a power-law in wavenumber, 
with spectral index n giving the scale dependence. Different 
inflation models predict different values for this spectral in- 
dex, and the high accuracy with which n can potentially be 
observed (Bond, Efstathiou & Tegmark 1997; Zaldarriaga, 
Spergel & Seljak 1997) promises v ery strong discrimination 
between different inflation models (Liddle & Lyth 1993; Lyth 



199C) 



This raises the question of whether future observations, 
in particular high-a.ccura,cv microwave background aniso- 
tropy measurements by satellites such as MAP and Planck, 
might be of such stunning accuracy that even the power- 
law ap proximation may prove inadequate. This has long 
been kjiown to be the case for 'designer' models of infla- 



just the point where perturbations on cosmologically observ- 
able scales are being created. Such models have long been 
regarded as rather unnatural, but retain interest as they 
are eminently testable by combinations of microwave back- 
ground and large-scale structure observations. 

A more pertinent question concerns whether the break- 
ing of the power-law might be detectable in the types of 
models regarded as theoretically most appealing. Such mod- 
els are of the slow-roll type (see Liddle & Lyth 1993 for a 
review). Several papers had pointed out that the spectral 
index would not be constant in typical models (Barrow & 
Liddle 1993; Copeland et al. 1993, 1994a), and an analysis of 
a rang e of slow- roll models was given by Kosowsky & Turner 
( 1995|) who wrote down a general formula for the scale de- 
pendence using the slow-roll expansion. Since their paper, 
however, theoretical prejudice has moved towards a subset 
of slow- roll inflation models, the so-called hybrid inflation 
mode ls ( |Linde 199l[ |Linde 1994 |Copeland et al. 1994b| ; |Lyth 



199f). These models have the disappointing feature obser- 
vationally that any gravitational wave contribution t o the 



micro wave anisotropies is expected to be negligible (Lyth 



tion (iKofman fc Linde 19871 ISalopck et al. 1989| ; [Hodgei 



et al. 199C; Hodges & Blumenthal 1990), where one tunes 



1997), but at least in some cases this appears to be compen- 



sharp features into the power spectrum by careful placement 
of strong features in the potential of the inflation field, at 



sated by a possible detectability of the scale dependence of 
the spectral index. 

An important question is the influence of the possible 
scale dependence on estimates of all other cosmological pa- 
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rameters, through the extra degeneracies introduced. We 
shall show that the cosmological parameters are only mildly 
affected. 



2 OBSERVATIONAL CAPABILITIES 

One of the crucial roles of microwave anisotropy experi- 
ments is to estimate cosmological parameters, ihese can 



be diviQed into two lypes. The cosmological parameters; 
such as the Hubble parameter, the density parameter and 
parameters describing the dark matter, give the evolution 
of the background space-time, the second set ol parame- 
ters describe the density perturbations which lead to the 
anisotropics. The simplest assumption, which coincides with 
the prediction of most inflation models and which we shall 
adopt throughout, is that the perturbations are adiabatic 
and Gaussian-distributed, in which case they are completely 
specified by the form of the density perturbation spectrum. 

It goes without saying that the more parameters we 
believe are needed to describe a data-set, the worse deter- 
mined they will be, so the input power spectrum plays a 
crucial role in the parameter estimation. Assuming at least 
that the power spectrum is suitably slowly varying, the best 
way to proceed (see e.g. Lidsey et al. 1997) is to expand its 
logarithm as a Taylor series in In k (where k is the comov- 
ing wavenumber) about some wavenumber ko, and truncate 
after some number of terms, i.e. 



InTMfc) 



lnVn(ka) + (n- l)ln 
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Here Vn(k) is the spectrum of th e cur vature perturba- 
tion, defined as in Liddle & Lyth ( 1993 ), which is scale- 
independent for a Harrison-Zel'dovich spectrum, which cor- 
responds to n = 1. Obviously ko is best chosen near the 
center of the range probed by the observations in question, 
but we will be more precise about this below. 

The Harrison-Zel'dovich spectrum amounts to truncat- 
ing this expansion after one term, and the power-law ap- 
proximation to truncating it after the second term. We shall 
investigate the possible inclusion of the third and fourth 
terms. 

Whether or not the power-law approximation is accu- 
rate clearly depends on the range of scales probed, measured 
by the size lnfc/fco reaches for a given data set. For Planck, 
which probes a range of multipoles from I — 2 to about 
2000, this term reaches about ±3 on either side of the cen- 
tral value. It is known that Planck can measure the spectral 
index n of a perfect power-law spectrum to an accuracy of 
better than ±0.01; however this does not directly lead to 
an estimate of the uncertainty in dn/dlnfc since this uncer- 
tainty also depends on the sensitivity of the given experi- 
ment to the power spectrum as a function of scale. Instead, 
we determine the errors on the first two derivatives of n by 
including them in the set of parameters to be determined 
from the data and using the Fisher matrix, as we now de- 
scribe. 

In principle, a brute-force maximum likelihood ap- 
proach would give the best estimates of a set of cosmolog- 
ical parameters from a CMB anisotropy data set, in the 



sense that the expected variances of the estimates are min- 
imized. In practice, such an approach would not be feasible 
for the quantity of data expected from the next generation 
of experiments (for example, satellite experiments will pro- 
duce CMB maps with at least 10 6 pixels). However, recent 
work on data analysis techniques suggests that near-optimal 
estimates of parameters coul d be achieved u sing computa- 
tionally fea sible algorithms (Tegmark 1997; Bond, Jaffe & 
Knox 1998). In any case, it is useful to consider what the 



smallest possible error bars on the parameter estimates will 
be and these can be calculated in a strai ghtforward way 



using the Fisher information matrix, atj (Tegmark et al 



1997). The covariance matrix for the parameter estimates is 
then the inverse of this matrix, and in particular t he stan- 
dard error for the estimate of parameter Sj is \/ (c~ 1 ^'-- 



This approach was used by Jungman et al. ( |l996a| ; |l996bj ) 
to determine the accuracy with which cosmological parame- 
ters could be determined from observations of temperature 
anisotropics, making use of a semi-analytical method to cal- 
culate the required angular power spectra and allowin g for a 
variation in the spectral index with scale. Bond et al. (1997) 
repeated their analysis, calculating the angular power spec- 
tra with a numerical Boltzmann code for greater accuracy, 
but did not consider deviations from a power-law pertur- 
bation spectrum. The improvement in parameter estimates 
when polarization inf orma tion is included was investigated 
by Zaldarriaga et al. (1997). We now generalize their results 
by including dn/dlnk and d 2 n/d(lnk) 2 in the set of param- 
eters to be determined. 

When both temperature and polarization anisotropy 
observatio ns are available the Fis her information matrix is 
given by ( ^aldarriaga et al. 1997 ) 



EE 



dC* 
dsi 



-Cov- 1 (Cxi,C Y i) 



dCy 

dsj 



(2) 



where Si are the parameters to be estimated and Cxi are 
the angular power spectra, with X and Y standing for T 
(temperature), E (even parity polarization), B (odd parity 
polarization) and C (E and T cross correlation) anisotropics. 
Cov -1 is the inverse of the covariance matrix between the 
estimators of the power spectra. The diagonal elements of 
this matrix are 
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and the non-zero off-diagonal elements are 
Cov(CtiCei) = , 2 1W C 2 ci 



Cov(CtiC C i) 



(21 + l)/ sky 
2 
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(4) 
(5) 
(6) 

(7) 



CcdCn + w^e 1 CT o) (8) 
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Cov(CeiCci) = ( 2 ; + 2 1)/k Cc^CEi+w^e 12 ^). (9) 
Here % and wp characterize the noise in temperature and 
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Table 1. Estimated parameter errors (one-sigma) for the Standard CDM model, under different 
assumptions for experimental configuration and generality of the underlying cosmological model. 



Parameter Planck HFI Planck 140 GHz channel 

no polarization with polarization 



5(C t )][ 2 /{C e )^ 2 


0.01 


0.01 


0.01 


0.005 


0.008 


0.01 


sn h h 2 /n h h 2 


0.006 


0.01 


0.01 


0.007 


0.008 


0.009 


Sn nl h 2 /h 2 


0.02 


0.02 


0.03 


0.02 


0.02 


0.02 


<5fi vac /i 2 /h 2 


0.06 


0.07 


0.08 


0.04 


0.04 


0.05 


<5t c 


0.12 


0.15 


0.15 


0.002 


0.002 


0.002 


6n 


0.006 


0.007 


0.007 


0.004 


0.004 


0.006 


Sr ts 


0.12 


0.14 


0.14 


0.04 


0.04 


0.04 


dn/dlnk 




0.02 


0.02 




0.009 


0.01 


d 2 n/d{lnk) 2 






0.02 






0.02 



polarization measurements respectively, / s ky is the fraction 
of the sky sampled and at is the beam width for the experi- 
ment. We calculate parameter uncertainties for two different 
specifications for the Planck satellite - one for four channels 
of the High Frequency Instrument (HFI) with no polariza- 
tion capability (using the experimental parameters of Bond 
et al. 1997) and one for the 140 GHz channel of the HFI with 
polarizers (using the parameters of Zaldarriaga et al. 1997) . 
In both cases we take the usable sky area as / s k y = 0.65. 

In addition to the specifications of the experiment be- 
ing considered, the estimated parameter errors also depend 
on the underlying cosmological model from which the data 
is drawn and on the choice of parameters to be determined 
from the data. For the latt er we adopt a subset of the vari- 
ables used by Bond et al. (1997) along with the derivatives 
of the scalar spectral index mentioned above. The densities 
of the various types of matter are specified by the parame- 
ters ujj = fljh 2 , where h is the Hubble parameter in units of 
100 km s" 1 Mpc -1 , and j=b, nr, vac refers to baryons, non- 
relativistic matter (baryons and cold dark matter combined) 
and the energy density associated with a cosmological con- 
stant, respectively. We assume there is no hot dark matter 
component and also that the Universe is flat, so the Hubble 
parameter is given by h 2 = cj nT +u> V!ic . The ionization history 
is characterized by the Compton optical depth tc from the 
redshift of reionization to the present, assuming complete 
reionization throughout that period. The normalization of 
the power spectrum is given in terms of the band-power 
correspondin g to the filter for the given experiment, {Ci}^ 2 
(Bond 1996). Finally, we consider the inflationary parame- 
ters: r ts which is the ratio of the tensor to scalar quadrupole 
moments, and the scalar spectral index n and its derivatives. 
They are discussed more fully in the next section. For the 
underlying cosmology we take a Standard Cold Dark Mat- 
ter (SCDM) model, with the following parameters; h = 0.5, 
fib = 0.05, fie™ = 0.95, fi vac = 0, T C = 0, r ts = 0, n = 1, 
dn/dlnk = d 2 n/d(lnfc) 2 = 0. 

The calculation of the Fisher matrix Eq. (Q) requires the 
derivatives of the angular power spectra with respect to each 
parameter, dCxi/dst. We determine these derivatives with 
finite differences, w ith the power spectra calculated using 
the CMBFAST code ( fceljak fe Zaldarriaga 1996 ; Zaldarriaga 



et al. IS 198). For some parameters these derivatives vary sig- 
nificantly within la of the target value, indicating that for 
these directions of parameter space the likelihood function 



is not well fit by a Gaussian. In these cases we use an itera- 
tive method such that the step size for each parameter used 
to calculate the derivative is approximately equal to its esti- 
mated la error. We have found that the resulting Gaussian 
approximation gives a good fit to the true likelihood func- 
tion, in particular correctly reproducing the strong correla- 
tions between some parameters. However, this dependence 
of the results on the details of the algorithm demonstrates 
that the given parameter errors should be taken as a guide 
only. Other uncertainties in the Fis her m atrix approach are 
considered by Magueijo & Hobson (1997). 

There is also freedom in the choice of the scale about 
which to expand the power spectrum, fco of Eq. (jl|). The es- 
timated errors on n and its derivatives are particularly sen- 
sitive to the choice of this scale, with values varying by an 
order of magnitude for different choices within the observed 
range of scales. However, an appropriate choice of fco can 
minimize these error estimates. In fact, when the true power 
spectrum is such that only n and dn/dlnk are non- vanishing 
it can be shownQ that fco can be chosen so that the error on 
n is unchanged from the case of a power-law spectrum, and 
is uncorrelated with dn/dlnk. This is because a scale exists 
at which the introduction of dn/dlnk does not contribute 
to the mean value of n, the mean being weighted by the er- 
ror estimation procedure. For the non-polarized configura- 
tion this optimal choice is fco = 0.064 Mpc -1 , while for the 
polarized configuration it is fco = 0.045 Mpc" 1 . Although 
in principle there is no guarantee that this is still the case 
when higher derivatives are present in the power spectrum, 
in practice we find that the above choices for fco give approx- 
imately the same errors on n and dn/dlnk when d 2 /d(lnfc) 2 
is included. We therefore give error estimates evaluated at 
these scales, bearing in mind that the errors at other scales 
can be much larger. 

The resulting parameter errors are displayed in Table |l|, 
for the two different experimental configurations and for 
increasingly general forms of the inflationary perturbation 
spectrum. For a power-law spectrum our estimates are in 
good agreement with those of other workers, and do not 
change substantially when the power-law assumption is re- 
laxed, except that as noted above the errors on n and its 
derivatives increase dramatically at scales away from fco. It 



We are grateful to Daniel Eisenstein for pointing this out. 
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should be remembered that these values are not exact one- 
sigma errors, given the uncertainties of the Fisher matrix 
formalism mentioned previously. 



3 SCALE-DEPENDENT SPECTRAL INDICES 
IN SLOW-ROLL INFLATION 

The results for the spectra from slow-roll inflation are by 
now extrem ely we ll known. We follow the definitions of Lid- 
dle & Lyth (1993). The spectrum V-n of the curvature per- 



turbation 1Z is given by 
'H\ 2 fH'~ 



8 V 



3mpj 



(10) 



Here V(cf>) is the potential of the scalar field driving inflation, 
and e is one of three slow-roll parameters defined by 



'"pi 
16tt 



m P1 V 
~8T ~V ' 

m% x V'V" 



64tt 2 V 2 



(11) 
(12) 
(13) 



where a prime indicates derivative with respect to the scalar 
field <j>. Note that e is positive by definition. Despite the 
square, £ 2 can be either positive or negative; the square is 
to indicate that it is second-order in the slow-roll expansion , 
which is an expansion in (m P id/d(j}) 2 ( |Liddle et al. 1994 - 
These parameters must be less than one for the slow-roll 
approximation to be valid. 

At a given point on the potential, corresponding to the 
location where the expansion in Eq. (Q) is carried out, the 
slow-roll parameters are in general independent since the 
derivatives can be freely chosen. There are further ones cor- 
responding to higher derivatives but we don't need them in 
this paper. 

The equations of motion for inflation give the slow-roll 
result 



d 



din k 



m% x YL A 

8?r V d(j> 



(14) 



which allows us to compute the spectral index and its deriva- 
tive with respect to the potential. The formulae are rem ark- 
ably simple; the spectral index is ( [Liddle fc Lyth 199$ 



n — 1 — 6e + 2r] , 



(15) 



(16) 



and its derivative is (Kosowsky & Turner 1995) 

_^- = -24e 2 + 16e f? -2e 2 . 
am k 

Along with these, another crucial inflationary observ- 
able is the influence of gravitational waves, relative to den- 
sity perturbations, on large- angle microw ave background 
anisotropies, given by (Liddle & Lyth 1992) 



rta 



14e. 



dn 



din k 



_ C 2 (grav) 
C2 (dens) 

Combining these equations gives 
: 0.12r t 2 s - 0.57r ts (l - n) - 2£ 2 



(17) 
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Figure 1. This shows contours of constant dn/dlnfc, assuming 
£ 2 = 0. Detection of the scale dependence at 95% confidence 
requires roughly | dn/d In k <; 0.01. 



which was first given by Kosowsky & Turner (1995). Re- 
member that in general £ 2 is independent of the other two 
parameters, and that it can be positive or negative. In Fig. 1 
we show the expected scale dependence under the assump- 
tion that £ 2 = 0; since in general it is unconnected to the 
first two this should be a good guide though obviously it is 
possible for some degree of cancellation to occur. For exam- 
ple, in the power-law inflation model the £ 2 term precisely 
cancels the scale dependence, though this is a very special 
situation. 

As we confirmed in Section 0, gravitational waves are 
only detectable (at 95% confidence level) if r ts ii 0.1 even 
i f cosmic microwave ba ckground polarization is measured 
( ^aldarriaga et al. 1997 ). The most optimistic scenario for 
measuring the scale dependence is if we can assume that all 
higher derivatives are negligible, and then with polarization 
a 95% confidence detection is possible if |dn/dlnfc| <; 0.01. 
Studying the first term of Eq. (n8|) , we immediately see that 
for it to be important then the tensors have to be detectable, 
by some margin. Much more interesting is the second term, 
because provided n is not too close to unity, it can give a 
detectable contribution to the scale dependence even if r is 
not itself detectable (though r ts cannot be very small). This 
is illustrated in Fig. 1; there exist regions where |dn/dlnfc| > 
0.01 but where r ts < 0.1. In fact, hybrid inflation models 
often have undetectable r ts but n significantly different from 
unity. The final £ 2 term in Eq. (^) is in general independent 
of the others and so nothing can be said directly, other than 
that a priori it is as likely to reinforce the scale dependence 
as partially cancel it. 

In conclusion, then, there exist significant parameter re- 
gions, including ones currently observationally viable, where 
the scale dependence of the spectral index can be measured. 
In such cases, it provides extra information concerning the 
inflation model which would not otherwise be available. In 
some cases, the scale dependence is measurable even when 
the tensor contribution is not. 
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4 SPECIFIC HYBRID INFLATION MODELS 

We end with a brief description of a particular type of in- 
flation model where dn/dlnk is detectable. Stewart (1997a; 
1997b) has recently presented two interesting models, the 
idea being to use quantum corrections to flatten the po- 
tential of the inflaton sufficiently to allow it to support an 
inflationary epoch. The argument is not difficult to repro- 
duce. Typical inflation models require a non-zero scalar po- 
tential energy V = Vo. In supergravity this induces soft 
supersymmetry-breaking masses for all scalar fields leading 
to a classical potential: 



V = V 



(19) 



with \A\ ~ 1. (Here the stringy Planck mass Mpi = 
mpi/\/87r has been set equal to one.) Unfortunately this 
does not lead to slow-roll inflation, since |V"/V| — \A\ ~ 1. 
Stewart pointed out that if <j> has either gauge or Yukawa 
couplings, to vector or chiral superfields with soft supersym- 
metry breaking masses squared of order Vo, then quantum 
corrections will renormalize the mass of <f> leading to a one- 
loop renormalization group effective potential of the form 



V = Vo 



l-i/(51n^ 2 + . 



(20) 



where It <tl 1 is the one-loop suppression factor, /(0) 



A + 0(8), with the expression being valid for Vo <C <t>" -C 1. 
Such a correction can lead to slow-roll inflation in certain 
regions of the potential. In particular, if <f)* is defined by 







(21) 



where /» = /(<51n 0*) and /»(x) = df,/dx then \V" /V\ ~ 
O(S) in the vicinity around 0* and the loop-corrected po- 
tential is flat enough to drive inflation even though the bare 
potential is not. The mechanism can work for f' r either pos- 
itive or negative. If it is less than zero then <^)* is a mini- 
mum of the potential and so a hybrid inflation mechanism 
is required to end inflation at some critical value </> c ; such 



a model is constructed in Stewart's second paper (Stewart 
1997b|!] I f fit > then cj>„ is a maximum of the potential. 
Stewart ( |1997e| ) considered such a model where </> rolled to- 
wards the true vacuum at <f> ~ 1. 

We are interested in the spectral index obtained and 
its slope. Introducing g(y) = —/*(») — S/2fl(x) where 
y = S\n((/> t /(j>) and defining 8 by S\n(j>* = —1, Stewart 
showed that, to lowest order in the slow-roll approximation, 
the spectral index n was given by 

n = 1 - 2g'J + 8(1 + A) ( 1 



exp 



-2 + 



y/l+ 3/4(1+ A) 



yj\ + 3/4(1 + A) 
- g'J(N - N b ) 



(22) 



where g', = 2(1 + A) /(A + y/A(A + l)), N - 7V fr is the 
number of e-folds of inflation from when the scale k oc e~ N 
leaves the horizon to (f> = <f>f r , the value of <j> when it begins 
to fast-roll down the potential and inflation comes to an end. 



Table 2. Values of n and dn/dlnk in different regions of param- 
eter space. 
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0.02 


20 


1.44 
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0.02 


40 


1.23 


0.008 
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0.04 


20 


1.18 


0.015 
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0.04 
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0.99 


0.005 



Since d/d\nk ~ -d/dN, it follows from Eq. (Q that 



dn 
dink 



= 8(1+^4) 1- 



giS exp 



a/1 + 3/4(1 + A) J 

g'J{N-N { , 



(23) 



y/l+ 3/4(1 +A) 



Typically the number of e-folds of inflation remaining 
after observable scales leave the horizon is 40 to 60, depend- 
ing on the energy scale and the reheating efficiency. However, 
this may be reduced by 10 or 20 if there is a second, low- 



energy period of inflation known as thermal inflation (Lytb 
& Stewart 1995). Low values of N — N{ r favour detectable 
scale dependence. 

The most natural scale for such a particle physics mo- 
tivated potential is Vo ~ M| USY where Msusy is the super- 
symmetry breaking scale in our vacuum. A value Msusy ~ 
10 10 -10 11 GeV corresponds to a gravity-mediated supersym- 
metry breaking. The corresponding COBE normalized value 
of (f>„ gives <f)„ ~ 10~ n , hence 8 ~ 0.04, a number consistent 
with that derived from a gauge coupling strength similar to 
that of the Grand Unified Theo ry gauge coupli ng inferred 
from LEP data, «gut ~ 0.04 ( Stewart 1997b ). However, 
there is significant flexibility in the allowed values of A and 
Vo, and hence S. A few values are tabulated in Table ^ to 
give a feel for the possibilities. What is encouraging is that 
provided there is a low number of e-foldings, the scale de- 
pendence can be detected by Planck or similar. 



5 CONCLUSIONS 

We have considered the likely magnitude and impact of scale 
dependence of the spectral index in inflationary cosmolo- 
gies. We have estimated the magnitude necessary to make 
the scale dependence detectable, and shown that slow-roll 
models of inflation, especially those of the hybrid inflation 
type, may give an observable effect. If so, this provides ex- 
tra information on the inflaton potential which would not 
otherwise be available. 

If scale dependence of n is considered, then there is a 
scale (around fco ~ 0.05 Mpc -1 ) at which the uncertainty 
in the determination of n is unchanged from the case of a 
power-law spectrum when the first derivative of n is intro- 
duced. We find negligible degradation of the error on n even 
when a second derivative is introduced, and have determined 
the anticipated errors on the first two derivatives of n at the 
preferred scale. We have found that at other scales the er- 
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rors on n are significantly greater; up to a factor of around 
10 greater for scales corresponding to the current horizon. 

The introduction of extra parameters to model the cos- 
mology implies that all the cosmological parameters will be 
more poorly determined. Fortunately, we have shown that 
the degradation of the uncertainty in these parameters is a 
small effect when the two derivatives are introduced. 



This paper has been produced using the Royal Astronomical 
Society/Blackwell Science DTfjX style file. 
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